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® A nuclear spectroscopy method for pulse height 
analysis of an electrical signal emitted by a radiation 
detector and including nuclear events, such as 
pulses, whose amplitude is a measure of the energy 
of the gamma rays collected by said radiation detec- 
tor, wherein (1) said signal is continuously converted 
to digital samples, at a given rate, and (2) each of 
the digital samples is processed so as to form a 
diaital image of each detected pulse. The energy of 



each pulse is calculated by summing all sample 
values representative of this pulse and the sample 
just preceding the first sample representative of a 
pulse, as well as the sample just following the last 
sample representative of the same pulse. 
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NUCLEAR SPECTROSCOPY METHOD AND APPARATUS FOR DIGITAL PULSE HEIGHT ANALYSIS 



BACKGROUND OF THE INVENTION 



1 . Field of the invention 

This invention relates generally to the analysis 
and process for indicating the amplitude distribu- 
tion of random-ampiltude pulse signals, such as 
those generated by a scintillation detector / 
photomultiplier. 
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2. Description of the prior art 

Such devices are commonly used in many 
technical areas where measurements involve nu- 
clear particles and radiation detection, one among 
others being e.g. the well logging techniques, 
wherein a tool is lowered in a well to carry out 
physical measurements. Of the many well logging 
instruments and techniques developed over the 
years to determine the hydrocarbon content and 
productivity of earth formations, the spectroscopy 
tool, by which energy spectra of the constituents of 
formation matrices and fluids are generated, has 
been found to provide information of particular val- 
ue in formation analysis. Typically, the energy 
spectra are obtained by detecting either natural 
gamma rays or gamma rays resulting from the 
interaction between formation nuclei and high en- 
ergy neutrons irradiating said formation, and con- 
verting each detected gamma ray into an electrical 
pulse whose amplitude is a measure of the gamma 
ray energy. These pulses are then sorted accord- 
ing to height in a pulse height analyzer to develop 
energy spectra characteristic of the constituents of 
the earth formations near the tool. 

The book "Radiation Detection and Measure- 
ment"< by Glenn F. Knoll (1979) depicts a typical 40 
pulse height analyzer, especially pages 720-725. 
The analyzer is disposed at the output of a scintil- 
lator detector linked to a photomultiplier, and is 
adapted to carry out an analog-to-digital conver- 
sion. Prior art analyzer is usually comprised of two 
main components, i.e. an analog/digital converter 
(hereafter referred to as ADC) and a memory. It 
also comprises, upstream of the ADC, an input 
gate which prevents pulses from reaching the ADC 
when the latter is busy, and another linear gate 
controlled by a single channel analyzer adapted to 
determine whether the amplitude is above a given 
threshold, and thus being representative of a pulse. 
Electronic filters, for amplifying and signal shaping 
nurnoses. are also provided at the input of the 
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pulse height analyzer. Since, basically, an ADC can 
only process one pulse at a time, simultaneous 
pulses, or successive pulses very close one to the 
other in time, cannot be processed, and thus de- 
tected; in other words, these pulses are lost, thus 
altering the reliability of the measurement. Accord- 
ingly, the faster the analyzer works, the more 
pulses are detected. However, the known analyzers 
comprise analogic circuitry, which implies (i) rela- 
tive low speed processing (ii) relative high costs 
and (iii) drift with time or temperature. The only 
attempt made until now for obviating the conse- 
quences of the low processing speed rely on dis- 
crimination circuits, which distinguish single pulses 
from stacked pulses, also called "pile-up" pulses; 
however, said circuits increase the complexity of 
the analyzer without giving a satisfactory answer to 
this problem. 

Furthermore, since such analyzers, depending 
upon the field of use, are located either downhole, 
in aircraft or in satellites, where room is limited and 
environmental conditions are severe, the pulse 
height analyzer should be compact and made of 
high reliability electronic components. Neverthe- 
less, known analyzers are made of numerous com- 
ponents, thus involving relative bulk, as well as 
difficult and costly adjustments due to functional 
dispersion. 

Moreover, a method has already been pro- 
posed wherein the analogic pulse signals are 
digitized and stored during a time window and 
wherein the digital data accumulated in the mem- 
ory are then processed by a computer. However, 
this known delayed process is time consuming and 
is not consistent with the requirements of the real 
time spectrum analysis technique. 



SUMMARY OF THE INVENTION 



A general object of the invention is to provide a 
pulse height analyzer showing high reliability, high 
speed processing, compactness, and being of low 
cost. 

The foregoing and other objects are attained in 
accordance with the invention by a nuclear spec- 
troscopy method for pulse height analysis of a 
signal, emitted by a radiation detector, such as a 
scintillator - photomultiplier device, and containing 
nuclear events, represented by pulses, whose am- 
plitude is a measure of the energy of the particles, 
such as gamma rays, collected by said radiation 
detector, said method comprising the following 
successive steps: 
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(1) continuously converting said signal to 
digital samples at a given rate: and 

(2) processing each of said digital samples 
so as to form a digital image of each detected 
nuclear event. 5 

The method of the invention also includes the 
step of detecting the arrival of a pulse, by compar- 
ing each incoming sample to a threshold value, so 
as to determine whether said sample is representa- 
tive of a pulse. More specifically, the difference w 
between the incoming sample value and the base 
signal which is free of nuclear events, is calculated 
and compared to said threshold. Furthermore, the 
time of arrival of any detected pulse is recorded. 

Moreover, said base signal value is continu- 75 
ously updated at each sample time arrival, so as to 
generate a current base signal value; preferably, 
the updated value is a weighted average of the 
incoming sample with the preceding sample, or 
preceding samples. 20 

Once a pulse has been detected, the energy of 
said pulse is calculated by summing the difference 
between each sample value (representative of said 
pulse) and the current base signal value, said sum 
being continuously stored in a register. 25 

During said energy calculation, the accumu- 
lated sum is compared to a preset value and said 
register is reset in case of overflow. 

Advantageously, for each detected pulse, the 
sample just preceding the first sample representa- 30 
tive of said pulse, as well as the sample just 
following the last sample representative of said 
pulse, are both taken into account for said pulse 
energy calculation. 

The method further includes the step of detect- 35 
ing a stack of successive pulses close one to the 
other in time. As a preferred embodiment, said 
detection step consists in making a count of the 
number of samples representative of a detected 
nuclear event and comparing said count to a pre- 40 
determined maximum count value. 

In another preferred embodiment, the sampling 
step is carried out by a flash analog-to-digital con- 
verter (ADC). Furthermore, the processing of the 
digital samples is carried out with a bit range 45 
greater than the bit range of the ADC, and math- 
ematical operations on samples (such as multipli- 
cation) are implemented through bit shifts. 

50 

BRIEF DESCRIPTION OF THE DRAWINGS 



FIG. 1 is a block diagram of an example of 
application of the invention in the form of a logging 55 
tool; 

FIG. 2 is a functional block diagram of the 
spectrum analyzer of the invention; 
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FIG. 3 shows a schematic representation of 
a pulse and the corresponding sampling times; 

FIG. 4 is a detailed flow chart of the pulse 
height analyzer of the invention; and 

FIG. 5 is a functional chart showing the de- 
tailed process of digital samples through the pulse 
height analyzer according to the invention. 

DESCRIPTION OF THE PREFERRED EMBODI- 

MENTS 



As mentioned, the present invention affords 
improved apparatus for pulse height analysis, 
which by way of example, has particular utility in 
well logging applications, and a representative em- 
bodiment of the invention is described and illus- 
trated herein in that environment. It will be under- 
stood, however, that this is merely illustrative of but 
one of the various applications for which the inven- 
tion is suitable. 

FIG. 1 illustrates a spectroscopy logging tool 
10 suspended in a borehole 12 by an armored 
cable 14. The tool 10 might be of the type de- 
scribed in U.S. patent 4,031 .367. issued to Richard 
D. Murphy on June 21, 1977. The cable 14 also 
provides electrical connections between the tool 10 
and the uphoie electronics, which include a com- 
puter 18, a visual recorder 20. a magnetic tape 
recorder 22 and additional intermediate signal pro- 
cessing circuits 16. Computer 18 suitably is a 
general purpose digital computer, such as that 
manufactured by Digital Equipment Corporation, 
Maynard, Massachusetts, under the designation 
PDP- 1 1 . The recorders 20 and 22 are coupled to 
the cable 14 by the usual mechanical linkage 24 
for making the customary records of the logged 
values as a function of tool depth. Appropriate 
power supplies (not shown) are provided at the 
surface and in the well tool, for supplying power to 
the uphoie and downhoie electronics. 

The tool 10 includes a pulsed neutron source 
26 which, under the control of a programmer 28, 
emits successive pulses or bursts of 14 Mev neu- 
trons into the formations surrounding the borehole 
12. Gamma radiation resulting from interactions 
between the neutrons and formation nuclei, e.g. 
inelastic collision gamma rays, is detected by a 
scintillation crystal 30. The usual shield 32 is inter- 
posed between the crystal 30 and the neutron 
source 26. The crystal 30 is optically coupled to a 
photomultipiier tube 34 which, in response to a 
detected gamma ray, generates a pulse signal 
whose amplitude is representative of the energy of 
the gamma ray. Output pulses from the photomul- 
tipiier 34 are passed by a spectrum analyzer 42 
including amplifying-filtering circuits 43 and a pulse 
height analyzer 44, the output of which is linked to 
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a memory 33, itself feeding a computer 35. 

At the surface, the data-bearing signals from 
the pulse height analyzer 44 are amplified, de- 
coded and otherwise processed in the signal pro- 
cessing circuits 1 6 for coupling- over a conductor 
bundle 54 to the computer 18. Digital spectroscopy 
outputs are transmitted to the tape recorder 22 
and, through appropriate digital-to-analog converter 
(DAC) circuits 58, to the visual recorder 20. Em- 
phasis is again laid upon the fact the present 
invention can also, as an alternative, be imple- 
mented in a natural gamma ray tool. 

FIG. 2 shows a diagram of the spectrum ana- 
lyzer 42 of the invention including amplifying-filter- 
ing circuits 43, receiving analog signals from the 
photomultipiier 34, and followed by the pulse 
height analyzer 44. Amplifying-filtering circuits are 
well known in the art and mainly comprise amplifi- 
ers and filters. Filtering allows one to shape the 
train of pulses in a predetermined fashion so that 
said pulses can be easily further processed. Filters 
are e.g. of the differential integrator type (CR-RC). 
and more specifically of CR-(RC) n type, realizing a 
Gaussian curve. More details about such filters can 
be found in the book (pages 61 4-623) by Glenn F. 
Knoll already referred to. 

FIG. 3 shows an example of a pulse 700 
(amplitude "A" versus time "t") having a Gaussian 
shape, while the essentially constant signal 710 
corresponds to the noise or background level, i.e. 
when no pulses are delivered by photomultipiier 
34. This constant signal is called the baseline B by 
one skilled in the art. Are also shown successive 

samples "a. b. c m, n, o n representative of the 

amplitude of the signal at successive instants ti - 
t,s 

Turning back to FIG. 2, the pulse height ana- 
lyzer 44 comprises an analog-to-digital converter 

60, which receives the analog signal issued from 
amplifying-filtering circuits 43 and converts said 
signal to digital samples at a given rate (ranging 
e.g. from a few MHz to a few tens of MHz). The 
ADC is preferably a flash ADC, such as, e.g., the 8- 
bit CMOS Flash converter HS 9584 manufactured 
by Hybrid Systems Corp. For each instantaneous 
sample S, the difference between said samples 
and the baseline B is compared to a threshold 
value T, arbitrarily chosen, via a threshold detector 

61. This step allows detection of an incoming pulse 
each time an incoming sample S, such as e.g. 
sample "e" of FIG. 3, of pulse 700 exceeds the 
baseline B by the threshold value T. 

in case the difference (S-B), i.e. sample minus 
baseline, is below threshold (T), said sample is not 
representative of a pulse, but of the baseline B. 
Said sample value- S is used for updating the 
baseline value B. via a baseline update unit 62, 
connected to the "NO pulse" output of said thresh- 



old detector 61, and the function of which will be 
further described. Any difference value (S-B) above 
the threshold T implies the instant sample S is 
representative of a pulse 700 , and is thus directed 

s from threshold detector 61 to a pulse processing 
unit 63. In most cases, the successive incoming 
samples, which are representative of a single pulse 
700, are directed via sample count detector 65 to 
an integration unit 64, which stores and processes 

70 said samples so as to calculate the energy of said 
pulse. 

However, since the kind of nuclear events the 
invention is directed to are random events, it often 
happens that successive pulses are quite simutta- 
75 neous or very close one to the other in time; thus, 
the corresponding nuclear event is actually a stack 
of successive pulses which involves loss of in- 
formation. This phenomenon is called "pile-up" by 
one skilled in the art. Means for detecting such 
20 pile-up events is dependent on the duration of the 
events, in other words, on the number of samples 
corresponding to the event. So as to detect and 
recover lost information, a given number of maxi- 
mum samples S for a single pulse, is chosen. In 
25 this case, said maximum number is set to seven. 
This maximum number is depending upon the du- 
ration of an average pulse and upon the sampling 
frequency. Accordingly, the pulse processing unit 
63 (at the output of threshold detector 61) is pro- 
30 vided wfth a sample count detector 65 (initials SC 
shown on FIG. 2 stands for Sample Count) which 
directs the incoming sample, depending upon the 
number of samples representative of the same 
pulse 700 which have already been processed be- 
35 fore said incoming sample, either to said integra- 
tion unit 64 (if SC is less than or equal to the 
maximum count 7) or to a pile-up unit 66 (if SC is 
more than maximum count 7). The output of pulse 
processing unit 63 and the output of said baseline 
40 update unit 62, are connected to a process next 
sample unit 67. The output of the process next 
sample unit 67 is fed to memory 33 and to the 
threshold detector 61. 

FIG. 4 shows a detailed flow chart of the pro- 
45 cess implemented in the pulse height analyzer 44 
according to the invention. For each new sample S 
issued from ADC 60. a processing value PV is 
calculated (block 401) as said new sample S minus 
baseline B. A pulse flag IF (block 402) is provided 
so and is SET each time an incoming pulse has been 
detected; said flag is otherwise NOT SET i.e. when 
no pulse has been detected. SET and NOT SET 
alternative is represented by two flow directions 
downstream of said IF flag, namely two threshold 
55 detection blocks 403 and 404. both performing a 
comparison between said processing value PV 
(sample S - baseline B) and the threshold T. Sup- 
posing the IF flag is NOT SET, threshold detection 
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403 leads to an alternative: if (1) the difference PV- 
T is negative, the baseline value B is thus updated 
(block 405) by weighted average of the incoming 
sample representative of the baseline and the pre- 
vious current baseline value, so as to generate a 
new current baseline value; and if (2) said dif- 
ference PV-T is positive, the incoming sample is 
thus representative of a pulse; the IF flag is SET 
(block 510), and the integration step 500 starts, 
while storing and gathering the different samples S 
representative of a pulse 700. The integration step 
500 allows also to take in account, as a pulse 
representative sample, the sample just preceding 
and the sample just following, respectively, the first 
and the last samples which are representative of a 
pulse. After the integration step 500 and the 
baseline update 400. the next sample is processed 
(block 600). Referring back to the IF flag 402. 
supposing it is in the SET position (i.e. a pulse has 
been detected), the incoming sample value S mi- 
nus the baseline value B. is compared to the 
threshold T (see threshold detection block 404) in 
the same way as already described, and a maxi- 
mum sample count detection is carried out (see 
blocks 406 and 407), so as to determine whether 
more or less than a given number of samples (in 
this case seven) have been identified as being 
representative of the same pulse. 

Supposing the threshold detection (block 404) 
is positive (i.e. the incoming sample is representa- 
tive of a pulse) therefore leading to the maximum 
sample count detector 406. A negative output of 
maximum sample count (from 406) leads to said 
integration 500; a positive output of maximum sam- 
ple count detection 406 means that the duration of 
the incoming pulse is longer than a given duration 
(corresponding to an average pulse) i.e., the incom- 
ing nuclear event is actually made of pile-up 
pulses. The incoming sample is thus directed to a 
"pile-up" processing unit 408, including a pile-up 
flag PF, which is moved to a SET position. The 
goal of said pile-up processing unit 408 is to re- 
cover the lost information from the different piled- 
up pulses, such units being well known in the art 
and examples can be found in the book by G.F. 
Knoll. The pile-up process is repeated for each 
sample fulfilling the corresponding condition. 

Supposing now the threshold detection 404 is 
negative, i.e. the incoming sample is not repre- 
sentative of a pulse, and more specifically, sup- 
pose that said incoming sample is the last of a 
sample series forming part of the pulse which has 
just been processed (since the integration flag is 
SET). Again, a maximum sample count detection 
(407) is carried out. A positive detection output 
from 407 means the pile-up process, which was 
ongoing, has to come to an end, and the integra- 
tion flag is thus moved to the NOT SET position. A 



negative maximum sample count detection from 
407 means the ongoing integration process is over 
and said negative pulse count leads to an addi- 
tional alternative, depending upon whether data 

5 overflow flag (OF) is either SET or NOT SET. If the 
OF flag is NOT SET. it is SET in the integration 
step 500 each time the pulse amplitude data 
reaches a maximum value corresponding to an 
overflow value. If the OF flag is SET, before pro- 

io cessing the next sample, the OF and IF flags are 
reset to the NOT SET position. It has to be noted 
that, at the very beginning of the pulse height 
analysis here above described, all flags are placed 
in the NOT SET position. 

75 Refer now to FIG. 5, which shows an example 

of implementation if the different functions, above 
described, performed by the pulse height analyzer 
44. 

The output of the ADC 60 is fed to a first adder 

20 71 which also receives the output of a complement 
circuit 70; said complement circuit delivers at its 
output the input signal with an inverted sign (e.g. 
positive value A thus becomes -A). Said first adder 
71 is linked to a shift register bank 72 which, also 

25 receives data from said complement circuit 70. A 
second adder 73 receives the output of said shift 
register bank 72 and of a first multiplexer 74, and 
feeds a second multiplexer 75, which directs data 
either to a first register 76 or to a second register 

30 77. First multiplexer 74 also receives data from a 
memory T which stores the threshold value T. The 
outputs of both first and second registers 76 and 
77 are fed to said first multiplexer 74. 

The way data are processed through pulse 

35 height analyzer 44, is schematically explained 
hereafter. As already mentioned, the update of the 
current baseline value B n *i . is carried out by 
weighted average, according to the following for- 
mula, given by way of example: 

40 B n +i = 1/4 S n ., + 3/4 B n [1] 

subscripts n n" and "n+1" standing for, respec- 
tively, previous and incoming values. Accordingly, 
the incoming sample Sn + 1 (which is representa- 
tive of the baseline) is granted a weight of 1/4, 

45 while the previous baseline value Bn is allotted a 
greater weight, of 3/4. Other weighting factors may 
of course be used, as well as the weighted average 
can also take into account more than one previous 
baseline value. 

so Multiplication of [1] by 16, and writing 12 as 

equal to 16-4, gives: 

16 B n ^ = (4 S n ^ - 4B n )+ 16 B n [2] 
where (4 S„* i - 4 B„) PV Le. 

the processing value already referred to in connec- 
55 tion with FIG. 4; 

Thus, divisions in equation [1] have been re- 
placed by multiplications. In the analyzer 44. mul- 
tiplications are carried out by bit shifts; e.g., the 
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product (4 S n +i) is made by adding to trie digital 
value Sn-i two binary zeros. The second register 
77 stores the value of the current baseline value 
(16 B n ), and delivers to said complement circuit 70 
the value (4 B n ). by removing two bits, in this case, 
the less significant ones. Said bits are used only 
for updating the baseline, and are not used during 
said integration step (corresponding to a sample 
above threshold). Second adder 73 performs the 
calculation according to formula [2]. The threshold 
detection, already referred to in connection with 
FIG. 4 (blocks 403 and 404), allows actuation of 
both multiplexers 74 and 75. Depending upon the 
result of the threshold detection, second adder 73 
is fed, by first multiplexer 74, with the value com- 
ing either from first register 76 or second register 
77. Also, depending- upon Jhejresult of said thresh- 
old detection, the value calculated by second adder 
73, is directed, through second multiplexer 75, ei- 
ther to the first register 76 for integration purpose 
(block 500 on FIG. 4), or to the second register 77 
for updating the baseline (block 405 on FIG. 4). 
Shift register bank 72 is designed to store the 
values preceding the incoming value (i.e. the value 
under process), so as to take into account for the 
calculation of the energy of the pulse under pro- 
cess, the sample just preceding the first sample 
representative of the detected pulse, as well as the 
sample just following the last sample representative 

of said pulse . 

While data issued by ADC 60 are digital data 
of 8 bit range, digital process in the pulse height 
analyzer 44 is carried out with a 10 bit accuracy in 
order to reduce round-off errors in the baseline 
subtraction. 

A final and very important benefit of the instant 
invention is that linearity error (both differential and 
integral) of the pulse height measurement, is sub- 
stantially reduced with respect to the linearity error 
of the ADC used. This is due to the fact that each 
pulse height determination is an accumulation of 
several ADC outputs, taken at different voltage lev- 
els in a random fashion, since there is no synchro- 
nism between the nuclear events and the clock 
used to sample the signal. Therefore, any local 
nonlinearities of the sample digitization process are 
reduced through averaging. 

Although the invention has been described and 
illustrated with reference to a specific embodiment 
thereof, it will be understood by those skilled in the 
art that various modifications and variations of that 
embodiment may be made without departing from 
the invention concepts disclosed. Accordingly, all 
such modifications are intended to be included 
within the spirit and scope of the appended claims. 

The digital circuitry of pulse height analyzer is 
preferably made as a custom integrated circuit. 



Claims 

1. A nuclear spectroscopy method for pulse 
height analysis of an electrical signal emitted by a 

5 radiation detector and containing nuclear events, 
represented by pulses, whose amplitude is a mea- 
sure of the energy of the particles, such as gamma 
rays, collected by said radiation detector, said 
method comprising the following successive steps: 

70 (1) continuously converting said signal to 

digital samples at a given rate; and 

(2) processing each of said digital samples 
so as to form a digital image of each detected 
nuclear event. 

75 2. The method according to claim 1 wherein 

the arrival of a pulse is detected by comparing 
each incoming sample value to a predetermined 
threshold value, so as to determine whether said 
sample is representative of a pulse. 

20 3. The method according to claim 2 wherein 

the difference between the incoming sample value 
and the value of the base signal, which is free of 
nuclear events, is calculated and compared to said 
threshold. 

25 4. The method according to claim 3 wherein 

said base signal value is continuously updated 
through each incoming sample, so as to generate 
an updated base signal value. 

5. The method according to claim wherein. 

30 once a pulse has been detected, the energy of said 
pulse is calculated by summing, the difference 
between each sample value representative of said 
pulse and said base signal value, said sum being 
stored as it is calculated. 

35 6. The method according to claim 5 wherein, 

for each detected pulse, the sample just preceding 
the first sample representative of said pulse, as 
well as the sample just following the last sample 
representative of said pulse, are both taken into 

40 account for said pulse energy calculation. 

7. The Method according to claim 2 including 
the step of determining whether an incoming de- 
tected nuclear event is a stack of successive 
pulses close one to the other in time, by making a 

45 count of the number of samples representative of 
said nuclear event, and comparing said count to a 
predetermined maximum count value. 

8. A nuclear spectroscopy logging method for 
determining characteristics of earth formations by 

so detecting particles (such as gamma rays) coming 
from the formations, resulting either from natural 
radioactivity, or from interactions between said for- 
mations with high energy neutrons irradiating said 
formations, and converting each detected particle 

55 into an electrical signal containing pulses whose 
amplitude is a measure of the energy of said 
particle, and by carrying out a pulse height analy- 
sis of said signal, said method comprising the 



11 EP0 396 



following successive steps: 

(1) continuously converting said signal to 
digital samples at a given rate; and 

(2) processing each of said digital samples 

so as to form a digital image of each detected 5 
pulse. 

9. Nuclear spectroscopy apparatus for pulse 
height analysis of an electrical signal emitted by a 
radiation detector and containing nuclear events, 
represented by pulses whose amplitude is a mea- io 
sure of the energy of the particles, such as gamma 
rays, collected by said radiation detector, said ap- 
paratus comprising means for continuously con- 
verting said signal to digital samples at a given 
rate, and means for processing said digital samples 75 
so as to form a digital image of each detected 
nuclear event through said digital samples. 

10. Apparatus according to claim 9 comprising 
pulse detection means able to compare each in- 
coming sample value to a predetermined threshold 20 
value and thus determining whether said sample is 
representative of a pulse. 

11. Apparatus according to claim 10 wherein 
said pulse detection means include means for cal- 
culating the difference between the incoming sam- 25 
pie value and the value of the base signal, said 
base signal being free of nuclear events, and for 
comparing said difference to said threshold. 

12. Apparatus according to claim 11 compris- 
ing means for continuously updating said base 30 
signal value through each incoming sample value, 

and for generating an updated base signal value. 

13. Apparatus according to claim 12 wherein 
said updating means includes means for making a 
weighted average of the incoming sample value as 
with the preceding sample value (or values). 

14. Apparatus according to claim 10 compris- 
ing means for calculating, once a pulse has been 
detected, the energy of said pulse by summing the 
difference between each said sample value which aq 
representative of said pulse and said base signal 
value, and means for storing said sum. 

15. Apparatus according to claim 14 wherein 
said calculating means includes means for compar- 
ing, during said energy calculation, the accumu- <s 
lated sum to a preset value, and means for reset- 
ting said storing means, in case of overflow. 

16. Apparatus according to claim 14 wherein 
said calculating means comprises means for in- 
cluding, for each detected pulse, the sample just so 
preceding the first sample representative of said 
pulse, and the sample just following the last sam- 
ple recresentative of said pulse, in said pulse en- 
ergy calculation. 

17. Apparatus according to claim 10 compris- 55 
ing means for determining whether an incoming 
detected nuclear event is a stack of successive 
pulses close one to the other in time, including 
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means for counting the number of samples re- 
cresentative of said nuclear event, and for compar- 
ing said count to a predetermined maximum count 
value. 

18. Apparatus according to claim 9 wherein 
said conversion means comprises an analog-to- 
digital converter (ADC) of the flash type. 

19. A nuclear spectroscopy logging apparatus 
for determining characteristics of earth formations 
by detecting particles (such as gamma rays) com- 
ing from the formations, resulting either from natu- 
ral radioactivity, or from interactions of said forma- 
tions with high energy neutrons irradiating said 
formations, and converting each detected particle 
into an electrical signal containing pulses whose 
amplitude is a measure of the energy of said 
particles, said apparatus comprising pulse height 
analysis means, including means for continuously 
converting said signal to digital samples, at a given 
rate, and means for processing each of said digital 
samples so as to form a digital image of each 
detected pulse. 

20. Apparatus according to claims 9 and 19 
wherein the digital circuitry of said pulse height 
analysis is made as a custom integrated circuit. 
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© Nuclear spectroscopy method and apparatus for digital pulse height analysis. 



LU 



.© A nuclear spectroscopy method for pulse height 
analysis of an electrical signal emitted by a radiation 
detector and including nuclear events, such as 
pulses, whose amplitude is a measure of the energy 
of the gamma rays collected by said radiation detec- 
tor, wherein (1) said signal is continuously converted 
to digital samples, at a given rate, and (2) each of 



the digital samples is processed so as to form a 
digital image of each detected pulse. The energy of 
each pulse is calculated by summing all sample 
values representative of this pulse and the sample 
just preceding the first sample representative of a 
pulse, as well as the sample just following the last 
sample representative of the same pulse. 
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